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Motivated by recent resistance data in high T c superconductors under magnetic fields paral- 
lel to the CuO layers, we address two key issues on the Josephson-vortex phase diagram, an 
appearance of structural transitions on the observed first order transition (FOT) curve in inter- 
mediate fields and of a lower critical point of the FOT line. The structural transitions reflected 
on the thermal transition curve in YBCO are identified based on our mean field calculation, and 
it is argued that a lower critical point more easily occurs in the present B 1 c case compared 
with that in B || c case. 
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§1. Introduction 

Josephson-coupled layered superconductors under magnetic fields parallel to the layers (i.e., in 
B 1 c) have at least two qualitatively different field ranges .!'& In low enough fields denned by 
p = 2ird 2 B\/ r T /(fio <C 1, the intrinsic pinning effect due to the layering is a weak perturbation and 
can be neglected for most purposes. Oppositely, in high enough fields p > 1, the layering effect 
is essential and strongly affects the p-dependence of fluctuation properties and of the transition or 
crossover lines. Here, d is the layer spacing, vT(> 1) the anisotropy of penetration depths, and 4>q 
the flux quantum. These two regimes are separated by an intermediate field range characterized by 
p ~ 1 in which, for instance, the transition or crossover temperatures defined in the disorder-free 
limit begin to lose their field dependences. Since this intermediate field range is reached by the 
applied fields of the order of tesla in BSCCO and underdoped YBCO, one needs to examine in 
details what are expected theoretically to occur in this field range of real systems including even 
some pinning disorder. Recent resistance data in high T c cuprates have suggested a first order 
behavior!' BSE! of the freezing (melting) transition in clean limit, the presence! Bi^ of its lower 
critical point (LCP) at least in YBCO (see below), and a reflection!' B of field-induced structural 
transformations between different vortex solids. 

In refs.2 and 8, we have argued that the freezing to a solid of Josephson vortices in disorder-free 
limit in B _L c is a first order transition (FOT) in all field ranges and that the FOT curve T m (B) in 
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high fields p > 1 becomes insensitive to p with increasing field. This is quite consistent with recent 
observations.! SB! We will focus below on the remaining two experimental findings mentioned 
above. First, to examine the reflection of structural transitions on the resistivity data, we need to 
clarify which types of vortex solids become ground states in clean limit in the relevant field range. 
As indicated in ref.2, in a narrow field range around a commensurate field where an intrinsically- 
pinned solid becomes the stable ground state due to a commensurability condition, the free energy 
in the solid phase is lowered reflecting a fluctuation-reduction due to the intrinsic pinning effect, 
and hence, the FOT temperature T m (B) is elevated. Since a unpinned solid or fractionally-pinned 
solid, i.e., a ground state with much weaker intrinsic pinning effect, can take the place of the pinned 
solid as the field slightly deviates from the commensurate field, an oscillating behavior of T m (B) 
may occur! and, in fact, seems to have been found in resistivity data of underdoped YBCO.Bi '! 
Although such an oscillating behavior was first speculated! for a putative vortex liquid-smectic 
transition, we do not believeii the presence of vortex smectic phase at least in p > 1 where T m 
should begin to become field-independent. Recent resistivity data in lower fields in BSCCo! have 
certainly shown a FOT between the vortex liquid and solid which, in contrast to the claim in 
ref.4, must not occur as the vortex smectic melting. This fact implies the absence of smectic phase 
even in p < 1. A key point is that the oscillating T m (B) curve was suggested in higher fields!!! 
p > 1 where the intrinsic pinning effect is essential. Such a high field FOT line with oscillating 
behavior is consistent not with the vortex smectic-liquid transition scenario! but with the liquid- 
solid transition scenario with field ranges,! in which the unpinned (floating) solid becomes the 
disorder-free ground state and hence, T m is lowered. According to the data,! only a couple of 
field windows of pinned solids are observed. Since, in high enough fields (p > 1.4), field-induced 
structural transitions no longer occur,! our theory may be able to clarify which vortex solids are 
realized in the observed oscillating melting behaviors. Calculated results on this issue will be given 
in §2. 

Another aspect relevant to the region p ~ 1, i.e., the presence of a LCP of the FOT curve, 
was suggested in resistivity data of YBCO with weaker fluctuation than of BSCCO: According 
to the underdoped YBCO datajli'! the in-plane resistivity {p a b) curves in lower fields suggest 
a continuous melting or glass transition, while p a b curves in higher fields have shown an abrupt 
decrease indicative of a FOT in the defect-free limit. Situation is similar to B || c case,0) in which 
a LCP has been detected in overdoped YBCO but not in BSCCO with stronger fluctuation and 
larger anisotropy. Therefore, it is natural to guess that the second order /5 a b-vanishing0) in tesla 
range of optimally-doped YBCO had been also a phenomenon in lower fields than LCP. Actually, 
no second order transition is suggestedBB* in the corresponding data in the low current limit, in 
both high and low fields, for clean BSCCO samples which have much stronger fluctuation and 
hence, a weaker pinning effect. Further, it has been quite recently clarified! that the out-of-plane 
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resistivity (p c ) in underdoped YBCO behaves in quite a similar manner to p a b and, for all fields, 
vanishes at essentially the same temperatures as those resulting from p a ^. Note that p c cannot 
vanish at a transition occuring without pinning disorder.il Therefore, the observed continuous 
vanishing!' of resistivity in lower fields is likely to be due to a Josephson- vortex-glass (JG) 
transition!!* occuring just above T m (B), while the almost discontinuous^!'!) resistive vanishing, seen 
in higher fields where the intrinsic pinning effect is essential, reflects a JG transition induced® by 
the FOT at T m (B). If this prediction is correct, an oscillating behavior suggested on both T m {B) 
curve and the glass transition line in high fields lying below T m {B) is not expected to appear on 
the continuous transition line in lower fields than LCP. However, a theoretical consideration is 
necessary in order to show that a situation in which the FOT is destroyed mainly by a LCP is not 
a rare event in real systems with random pinning effect under B 1 c and will be performed in §3. 
Finally, in §4, our theoretical results are compared with available resistivity data in B _L c. 

§2. Sequence of Structural Transitions in p ~ 1 

In this section, we reexamine the field-induced sequence of structural transitions between Joseph- 
son vortex solids, particularly in the intermediate fields around p ~ 1. Structural transitions in 
this region were considered only briefly in the previous works& (see Fig. 7 (a) in ref.2) and will be 
compared in details on this occasion with available data.BII* In addition, a recent study© using the 
London model with uniaxially periodic potential has suggeted an importance of rotated (pinned) 
solids in B _L c in which all vortices are pinned between two neighboring layers. Here, we primarily 
examine energies and stabilities of rotated solids possible around p ~ 1. In §4, it will be pointed 
out that a couple of rotated solids are good candidates of ground states reflected in the observed 
oscillating melting (and glass) curves. 

Following ref.2, we apply the LLL approximation of order parameter ip, defined!!* in continuous 
x-y space, to the Lawrence-Doniach model 



Hld = J dy J dx ^ exp | i 27T ™ x j 



e \iP(x, y)\ 2 + $\(-id y + r B 2 x)i/j(x, y) )| : 



+r- 1 d-%V(^, V) -^(x + d,y)\ 2 + h -\^{x, y)\ A ( 2. 1 ) 

rewritten!) in terms of the Poisson-summation formula, where £o is the in-plane coherence length, 
rs = \J 4>o/(2ttB), b > a constant, and Eq ~ (T — T c q)/T c q. Spatial variations parallel to B were 
neglected by focusing on the mean field solutions. Further, in the x-y plane perpendicular to B, the 
out-of-layer direction and the in-layer direction are denoted as x and y-directions, respectively, and 
the type II limit with no variation of flux density was assumed. In this representation, the "core" 
of a Josephson vortex is expressed as that of the ordinary singular vortex. Since the thickness of 
each superconducting layer is infinitesimal in the Lawrence-Doniach model, all vortex "cores" are 
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situated between the superconducting layers. In a single vortex picture, each vortex is expected to 
sit in the middle of an inter-layer spacing. However, the type II limit assumed here cannot reduce 
to such a single vortex picture because the interaction-range between the field-induced vortices is 
infinite in type II limit, and hence, the vortices in type II limit may form a solid structure favored 
by the long-ranged interaction between them without sitting in the middle of inter-layer spacings. 
Actually, it is found!) that fractional pinned solids, in which a fraction of vortices are not located 
in the middle of inter-layer spacings, can become ground states in clean limit at fields determined 
by their commensur ability conditions. 

A periodic solution including rotated lattices is expressed, within LLL in this representation, in 
the form 

V>(>> V) = V rl/4 — i7? Yl exp(i/cmy - 7- n-(x + kmr 2 B ) 2 + — to 2 ), (2.2) 

V 7T ' m 2r B R 

where R is a rational number. Although several rotated (pinned) solids in which R takes integer 
values larger than 2 were studied in ref.2 (see the sentences below eq.(2.4) there), they were not 
identified with ground states, possibly because of too low field values. As shown below, rotated 
structures with a fractional R can have energy remarkably lowered due to the vortex confinement 
induced by the layering and hence, become the ground states in the intermediate fields p ~ 1. 
Although the rotated version of the fractional pinned solids can also be created by choosing a 
fractional R, most of them seem to be fragile in real systems with disorder^ (see §4). Here we 
will focus on strongly (intrinsically-) pinned solids called in ref.2 integer pinned solids, in which 
each vortex is located in the middle of an inter-layer spacing. An integer pinned solid is defined 
by kr\ = wd with a positive integer w and has (w — l)d vacant inter-layer spacings between two 
neighboring occupied inter-layer spacings. Namely, in a w = 1 solid all inter-layer spacings are 
occupied by vortices. Further, we will not consider here integer solids with w > 2 but will focus on 
rotated w = 1 solids for a reason to be clarified in §4. 

Based on eq.(2.2), the rotated w = 1 solids are created as follows. The nonrotated w = 1 
solid is represented in terms of the parameters, R = 2 and kr B = d, and have the basis vectors 
= (27r/k)y and a^ 2 ) = (ir/k)y — wdx (see Fig.l (a)). Now, let us rotate these basis vectors 
counter-clockwise by an ang le 6, a® -> ak')(0). It is easy to find that the condition under which 
the resulting rotated solid is intrinsically pinned by the layers is Mok\o) = Na,x (0) = MNd, 
where M and ./V are positive integers. The resulting structure is represented within LLL in terms 
of eq.(2.2) by setting kr 2 B = d and 

2M-N K ' 

and will be called below (M, N) state. The commensurability condition of this structure is 

P rl M 2 + N2-MN' 1 ' 
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which means that the state (N, M) is denenerate in energy with the (M, N) state. Actually, 
it is seen that they can be transformed to each other by reversing the field direction B —* — B. 
According to the expressions mentioned above, the rotated states with JV = 1 or M = 1 belong 
to w = 1 solids with no vacant inter-layer spacings. Since a more energy gain due to the intrinsic 
pinning is expected to be maximized in a w = 1 state, we will primarily consider states with TV or 
M = 1. Actually, for instance, the (2, 3) state is equivalent to the (3, 1) one. Further, as shown 
previously,!^ the commensurability conditions for nonrotated integer solids with (w — l)d vacant 
spacings in the x-direction are p = ^ir/w 2 and ir/(V3w 2 ). From this and eq.(2.4), the (M, M) 
and (2M, M) states are found to be merely alternative representations of nonrotated solids. 

The rotated w = 1 solid to be realized in a highest field is the (1, 3) state, illustrated in Fig.l (b), 
or its reversal one (3, 1) (Fig.l (c)), and their commensurate field is given by p\% = \^3ir/7 ~ 0.77. 
The rotated state realizable in the next lowest field, with decreasing field, is the (1, 4) or (4, 1) 
state which becomes commensurate when pu = \/3n/13 ~ 0.42. Since the layering is not effective 
with decreasing field, other rotated states occuring by lowering field further will not be considered. 

Keeping these commensurate fields for various structures in mind, the sequence of ground states 
changing with increasing p can be determined!) by computing the value of Abrikosov factor fix 
of each structure for various p-values. The resulting sequence of ground states is given by the 
table of Fig. 2, where we take account only of the w = integer nonrotated solids and (M,l) rotated 
solids as candidates. In the p-ranges shown in blank spaces of the upper row of Fig. 2, fractional 
pinned solids or the unpinned (floating) solid become ground states. It is found that the w > 1 
nonrotated integer solids become ground states just at their commensurate field values, while a 
(M, 1) state, particularly the (3, 1) state, is the lowest in energy over a broad range including 
its commensurate field. In the case of (3, 1), /3a < 1-15959 (the value for the floating solid) in 
the range 0.725 < p < 0.96 including its commensurate field p = 0.77, and j3\ reaches its lowest 
value 1.157 near p = 0.8. Further, we found that the (4, 1) state becomes the ground state in 
0.416 < p < 0.421. This result that, at least in intermediate fields specified by p-values of order 
unity, a rotated w = 1 solid is minimized over a broad range including its commensurate field value 
is a reflection of the intrinsic pinning effect enhanced by the absence of vacant inter-layer spacings. 
Thus, it may be possible that, if the model (2.1) is analyzed beyond the LLL approximation, a 
rotated w = 1 solid becomes the ground state at a commensurate field of a w > 1 nonrotated solid. 
Actually, we have found in ref.2 that the w = 1 nonrotated solid is lower in energy than the w = 2 
one when p = \/3vr/2 2 . 

It should be noted that the period in the x direction of the (M, 1) state is 2(1 + M 2 — M)d, 
while the corresponding period of the w = 2 nonrotated solid is Ad. Hence, to find the (M, 1) solid 
in a numerical simulation in fields near its commensurate field, setting the system size in the x 
direction to a multiple of 2(1 + M 2 — M)d is necessary. Otherwise, a numerical simulatioin would 
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create a wrong ground state at lower temperatures. 

Before closing this section, it is necessary to clarify whether a superposition of the two states, 
(M, 1) and (1, M), results in a lowering of energy, since they are degenerate in energy with each 
other. In fact, the authors of ref.13 have argued that a mixing of the two pure states will be a 
representation of a waving vortex solid structure. On the other hand, the present author has given 
an evidence of such a waving solid in a limited field range (1.26 < p < 1.4) by showing^ that this 
structure is created due to a shear instability of the w = 1 nonrotated solid. To examine whether 
the former mechanism© is possible, we have computed the energies of various superpositions of 
the (1,3) and (3, 1) states in the range 0.72 < p < 0.93 (see Fig. 2) and found that such a mixing 
always raises the energy. Namely, in fields around the commensurate field, either the (1,M) or (M, 
1) states is spontaneously realized as the ground state. Actually, we have checked that the shear 
modulus of the (1, 3) state is positive in the above-mentioned field range, guaranteeing the stability 
of pure rotated states. It is physically natural that a mixing of the two pure rotated states does 
not result in a waving ground state solid: The picture on the waving solid argued in ref.13 requires 
that, when moving along the x direction (_L layers), the (1, M) and (M, 1) states alternate with 
domain walls intervening between them. Due to such domain walls with positive energy, a pure 
rotated solid should have a lower energy compared with any mixing of two equivalent pure rotated 
solids. 

§3. Presence of a Lower Critical Point 

In ref.8, we have studied the JG transition in B 1 c in higher fields (p > 1), where the FOT 
curve T m in disorder-free case is insensitive to p, and the glass transition temperature Tjg is lower 
than T m (p). However, a remarkable deviation of Tjq(p) from T m (p) to lower temperatures has been 
seen only in dirtier BSCCO under higher fields^) where the FOT in B 1 c was not detected. In 
underdoped YBCO with much smaller anisotropy than in BSCCO, the TjG(p)-line in p 3> 1 does 
not deviate remarkably from T m (p) particularly in cleaner samples^ (see also §4). 

Rather, a continuous vanishing of in-plane resistivity p a b in YBCO has been often observed in 
lower fields where T m {B) is sensitive to B. Taking account of the corresponding p c -vanishing^ this 
continuous transition is, just like that explained in B || c case, likely to be the JG transition!) 
occurring below a LCP. Further, since no discontinuous resistive vanishing occurs in the same 
sample!) in B || c, we feel that the FOT in B 1 c case tends to be destroyed, in contrast to B || c 
case, primarily from the lower field side by creating a LCP. To show that this is indeed so, we 
consider an extrapolation to higher fields of a JG curve computed in p <C 1 which, just as the 
thin solid curve shows in Fig. 3, lies in p < p* above T m (B) due to point disorder with moderate 
strength. Based on the familiai@ results in B || c case, one might expect the Tjq(B) line to have 
a functional form similar to T m {B) (described as the chain curve in Fig. 3) even in B _L c case 
and hence, to become insensitive to B in higher fields. If so, the FOT would not occur at any 
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field because the resulting JG transition curve would always lie above T m (B), and the fact, that 
Tjg > T m in lower fields, would imply that the FOT be destroyed mainly through a reduction of 
an upper critical point. As is seen below, however, this picture does not apply to the B 1 c case. 

Upon cooling from higher temperatures, the JG transition is signaled by a divergence of the 
glass susceptibilityEl'li' xg calculated in terms of the model (2.1), and the diagrammatic form of 
Xg is typically a series of ladder diagrams composed of bare disorder (random pinning) lines. As 
shown in B || c casejll* the vertex correction to each bare disorder line arising from a precursor of 
vortex solidification enhances the disorder strength and Tjg , while the vertex correction originating 
from a high temperature fluctuation effect unrelated to the vortex solidification reduces Tjg- The 
former effect in high field limit for B 1 c case was studied in ref.8. Here, we consider, by assuming 
Tjg > T m in lower fields, just the latter effect which can be represented by the familiar diagrams of 
RPA type. In low fields where the layer structure is negligible, the XG- ex P ress i° n is essentially the 
same as that for the isotropic 3D GL model and, in LLL approximation, gives a TjG-curve yielding 
the LLL scaling© and may lie above T m (B) for a moderate disorder strength. Since it is difficult 
to directly perform the corresponding calculation in the intermediate fields, let us consider here the 
high field limit (p > 1), in contrast to in ref.8, by neglecting the vertex correction related to the 
vortex solidification effect. As is intuitively clear, the model (2.1) loses all field dependences in this 
limit and takes the same form as the GL model in 2D and zero field except for the shift T c q—T C 2 of the 
position of its mean field transition (see §4). On the other hand, it is well knowr© through formally 
the same calculation as that perfomed in the context of spin-g lass0> that a superconducting- 
glass transition in homogeneous (i.e., nongranular) materials under zero field does not occur even 
at the mean field level (i.e., irrespective of dimensionality). Actually, since, in the irreducible 
vertex in the ladder representing xg, the fluctuation vertex correction to the disorder strength 
outweighs the temperature dependence accompanying the bare disorder strength, the irreducible 
vertex itself decreases on cooling and consequently, xg cannot diverge at any temperature even 
if the renormalization of glass fluctuations, leading to a dimensionality dependence, is neglected. 
This result Tjg(p — > oo) — > obtained without effects of vortex solidification, when combined with 
the result on Tjg in ref.8 described in Fig. 3 as the thick solid curve, implies that, as described by 
the solid curves in Fig. 3, any TjG-line lying above T m in lower fields must cross T m (£?)-curve at 
some field and lie at lower temperatures than the disorder-free T m (i?)-curve (the chain curve in 
Fig. 3) in higher fields. Clearly, a point essential to reach this conclusion is that, in contrast to in 
B || c case, the T m (i?)-curve becomes insensitive to p in p > 1, i.e., that T m drastically changes its 
p-dependence in an intermediate field p ~ 1, while Tjg has no such a property. Therefore, a LCP 
inBlc case is expected to, in many cases, be created in the range p ~ 1 so that the FOT in p > 1 
with strong intrinsic pinning effect remains. 
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§4. Discussion 

Based on our results in the preceding sections, let us discuss here the doping dependences of 
resistivity data in ref.4 and the resulting phase diagram in B 1 c. As mentioned in §1 (see also 
the comment!' on ref.4), the oscillating behavior on the vertical portion of the FOT line is a clear 
evidence that the transition corresponds not to the putative liquid to smectic transition but to 
a liquid to solid (strictly speaking, liquid to slushlllHl^ ) transition argued in ref.2. In high fields 
p>l, this liquid to solid (slush) transition line T m (p) was denoted in ref.8 as T sc and, according 
to §5 in ref.2, is given by 

T m {p » 1) * T c0 (l - ^r- 1 )^ + Coo eg))" 1 (4.1) 

(see the paragraph including eq.(5.1) in ref.2), where e^q = 167r 2 k-BT c o\ 2 (0)/((pQd) is 2D Ginzburg 
number, A(0) the in-plane penetration depth, Coo(> 0) is a constant of order unity, and the com- 
bination T C 2 = T c o(l — 2r _1 (^o/^) 2 ) in eq.(4.1) expresses the mean field transition temperature in 
p —* oo limit.!' Now, let us compare eq.(4.1) with Fig. 4 in ref.4, which shows that the transition 
temperature vertical (i.e., independent of B) in higher fields increases with underdoping. Since 
both r and A(0) increase with underdoping, while a doping dependence!) of T c q is negligible in the 
doping range of samples in ref.4, the doping dependence of the vertical transition line is primar- 

(2) 

ily due to the T-dependence of eq.(4.1), i.e., of Namely, the doping dependence of e G may 
be neglected compared with that of T. We emphasize that, in the phase-only approximation of 
the Lawrence-Doniach model, T C 2 is obtained as T c q, and hence that the above-mentioned doping 
dependence of T m (jp > 1) cannot be explained in terms of phase-only mo delsfl) such as the XY 
model. 

Next, let us consider at what field T m {B) begins to become insensitive to B. To clarify this, 
it will be sufficient to see at what field, called I?x hereafter, the portions of T m {B) in higher 
{p S> 1) and lower {p <C 1) fields meet when they are interpolated with each other. Since eq.(4.1) 
is independent!' of an approximation method, while the portion in p <C 1 where the layering is a 
secondary effect depends highly on the approximation methods, we need to consider under what 

condition the LLL approximation for the p < 1 region becomes inapplicable. Estimating a crossover 

(3) 

field Bcr between the LLL and the London behaviors in B || c case was performed in ref.16 (see 
eq.(4.1) in ref.16). Within the anisotropic 3D GL model,© a crossover field in B 1 c corresponding 
to this is easily found, and the LLL approach is semi-quantitatively valid when 

B > VTb£\ (4.2) 

(3) 

In ref.16, B£ J was estimated to be 1.7 tesla for the optimally-doped YBCO. Using this and based 
on eq.(4.2), the r.h.s. of eq.(4.2) may be, due to its T-dependence, beyond 50 (T) for underdoped 
YBCO even if the doping dependence of A(0) is neglected. Through this fact, the T m (i?)-curves in 
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p < 1 for underdoped YBCO will be assumed to be well described in the London limit of anisotropic 
3D GL model. The corresponding T m (B) is nothing but T^ip) curve given in ref.8 and has the 
form Tm — T C 2(p)(l — Eq y/p). By comparing this with eq.(4.1), Bx appropriate to underdoped 
YBCO is defined assuming p to be a constant, i.e., Bx oc T -1 / 2 , where d was assumed to be 
doping- independent. By analyzing Fig. 4 in ref.4 in terms of this Bx and eq.(4.1), we find that, for 
each sample in ref.4, y/T — 32 (T750), 14 (T700), and 8 (T650), by assuming other parameters 

(2) 

to be doping-independent. Consistently with these T values, we find CqqEq ~ 0.043. It means 
that, when taking a reasonable value A(0) ~ 10 3 (A), the number Cqo is about 2.0 which coincides 
with the numbei@) giving the Kosterlitz-Thouless transition point in 2D films in zero field. This 
is reasonable because the present problem would reduce to the 2D system in zero field if the shift 
T c q — T C 2 of the mean field boundary be absent. 

We next try to describe the oscillating T m (B) of the sample T750 in ref.4 by ascribing its origins 
to a couple of pinned ground states. Figure 2 in ref.4 shows that the oscillating T m and Tjq lines 
suggested from the resistivity data include three stability regions of pinned solid states in fields 
below 14 (T). The FOT at T m {B) in clean limit seems to be reflected in the data above 4 (T) 
and below 12 (T). Namely, both upper and lower critical points should be present there: We note 
that the resistive vanishing at 14(T) in Fig.l of ref.4 is clearly continuous and that the presence 
of a upper critical point below 14 (T) is expected. Of course, the feature that the transition 
curve is vertical (5-independent) above 4 (T) originates from the intrinsic pinning effect due to 
the layering. One of remarkable features is that the middle window, suggestive of a pinned ground 
state structure, around 6 (T) is remarkably broad. First, the large T (~ 1 x 10 3 ) value of the sample 
mentioned above implies that the p ~ 1-value for this sample corresponds to about 10 (T). Since, 
according to the present Fig. 2, no structural transition will occur in p > 1.26 of real systems ,0 the 
suggested field-induced transition near 10 (T) will imply an onset of the w = 1 waving solid lying 
near p = 1.26. Then, turning to lower fields, it is quantitatively reasonable to assume that the 
wider window, suggestive of a stronger intrinsic pinning, around B ~ 6(T) will correspond to the 
window in Fig. 2 in which the (1,3) or (3, 1) state is stable. As mentioned in §2, the rotated w = 1 
solids such as a (M, 1) state are favored, even if the field deviates from their commensurate field, 
by the confinement effects of the layering and may have their wide widths in which these solids have 
the lowest energy. This tendency is lost with increasing M. Hence, a width of a stability region of 
a (M, 1) state relevant to lower fields should be narrower. If the w > 1 solids are not realized, the 
third pinned region near 3.3(T) suggested in Fig.l of ref.4 will be a reflection of the (1,4) (or (4,1) 
) state. Although the commensurate field p = 0.42 of the (4,1) state suggests that, in the data of 
ref.4, this state should be realized in slightly lower fields, and the w = 2 nonrotated solid is more 
reasonable in view only of commensurate field values, a rotated solid may be minimized, due to the 
vortex confinement effect, at a commensurate field of a w > 1 nonrotated solid (see §2). In fact, 
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although aw = 3 nonrotated solid has its commensurate field p = 0.59 and hence, is expected to be 
realized near 5(T) in Fig. 2 of ref.4, there is no sign of realization of the uu = 3 state in data of ref.4, 
possibly reflecting its disappearance!) due to the random disorder. Through this fact, there will be 
no reason why the w = 2 state should be realized in the same real sample. Therefore, we expect 
that the third pinned state near 3.3 (T) will correspond to the (4,1) (or (1,4)) state belonging to 
the w = 1 states. 

Here, we comment on the positions of Tjq(B) above an upper critical point (UCP) of underdoped 
YBCO according to the result in ref.8. Broadly speaking, Tjq(B) above UCP is less steep with 
increasing B, and the JG-transition curve derived in ref.8 

T*(A) - T 

B JG (T)^B ^. , (4.3) 

becomes applicable with increasing B, where Bq = 0o/(27rd 2 vT~ £q ), and A denotes the random 

(2) 

pinning strength. Since, assuming that only T is doping-dependent and that e G — 0.02 and 
d ~ 10(A), B for YBCO is estimated to be W 4 /VT tesla, the B JG (T) curve for underdoped 
YBCO in ref.4 with T < 10 3 is found to be more steep than the corresponding H C 2(T)-lme in B || c 
case. This is in contrast to that in BSCCO which was primarily considered in ref.8 and found to 
be much flatter in the B-T diagram. 

It will be necessary to explain why an oscillating behavior of T m (B) and its LCP were not 
observed in the corresponding data!'@) of BSCCO. The absence of LCP in BSCCO under B _L c 
does not contradict its presence in YBCO. As in B || c case its origin can naturally be ascribed 
to a stronger fluctuation and a much larger anisotropy in BSCCO both of which reduces the LCP. 
Further, since the FOT in B || c is present in the sample in ref.7 in contrast to the YBCO sample 
in ref.5, a weaker disorder is anticipated for the sample,^ which also results in a reduction of LCP. 
Rather, it should be noted that structural transitions of Josephson vortices have not been reflected 
at least in the vertical portionB of the resulting FOT curve. Instead, a couple of structural transition 
fields suggested from the datai in tilted fields were found to remarkably decrease with increasing 
temperature. By examining the BSCCO dataHEt altogether, the structural transition fields near 
T m (B) likely lie below Bx (see a preceding paragraph), i.e., in lower fields than the vertical portion 
of T m (B). This observation cannot be explained within the mean field approximation and suggests 
a strong T-dependence of the structure transition fields. This strong T-dependence may be due to 
an effective reduction of anisotropy occuring from the fact that the thermal fluctuation in BSCCO 
is not negligible even below T m (B). To resolve this issue, a quantitative study of fluctuation effects 
in the solid phases is needed and left to a future work. @ 
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Figure Captions 
Fig.l 

(a) Basis vectors a^ 1 ) and a^ 2 ) in the w = 2 nonrotated solid described in the x-y plane, where 
the horizontal lines denote the superconducting layers, and the open circles indicate the "cores" of 
Josephson vortices, (b) Vortex configurations in the w = 1 rotated solid (3,1), and (c) the one in 
the (1,3) state. 



Fig.2 

Table expressing the sequence of ground states at low enough temperatures. Here, the states 
denoted as w = integer imply the nonrotated integer pinned solids, and (M,l) states are rotated 
w = 1 solids. The blank regions on the upper row, such as the one including the p = 1 value, imply 
the field ranges in which the unpinned or fractionally pinned solids become ground states. Note 
that the structures expected to occur in p < 0.4 are not shown. 

Fig.3 

Schematic phase diagram in B _L c. Thin solid curve implies the Tjq(B) curve obtained by 
assuming Tjq(B) > T m (B) (the chain curve), and hence, only its portion in p < p* (i.e., below 
the LCP) plays the role of true Tjq(B). The true Tjq(B) curve in p > p* , examined in ref.8, 
corresponds to the thick solid curve. The thermal FOT cannot© occur on the portion in p < p* 
of the chain curve. 
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